Nature of fault planes in solid neutron star matter 
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ABSTRACT 

The properties of tectonic earthquake sources are compared with those deduced here for fault 
planes in solid neutron-star matter. The conclusion that neutron-star matter cannot exhibit 
brittle fracture at any temperature or magnetic field strength is significant for current theories 
of pulsar glitches, and of the anomalous X-ray pulsars and soft-gamma repeaters. 



Subject headings: stars:neutron - pulsars:general - X-rays:stars 



1. Introduction 

It is widely assumed that brittle fractures 
caused by Maxwell or other stresses in neutron 
star crusts are involved in a number of phenomena, 
for example, the soft gamma repeaters (SGR), the 
persistent emission of the anomalous X-ray pulsars 
(AXP) (Thompson et al 2000), and large pulsar 
glitches (Ruderman, Zhu & Chen 1998). In cur- 
rent theories of the AXP and SGR sources, brittle 
fractures, propagating with a velocity of the or- 
der of the shear- wave velocity c s , generate shear 
waves which in turn couple with magnetosphcric 
Alfven modes. At angular frequencies u> ~ 10 4-5 
rad s _1 , the coupling is thought to be an efficient 
mechanism for energy transfer to the magneto- 
sphere, as shown by Blaes et al (1989). Statistical 
comparisons of SGR burst properties with those 
of terrestrial earthquakes are not inconsistent with 
the brittle fracture assumption (Hurley et al 1994, 
Cheng et al 1996, Gogus et al 1999). However, the 
purpose of this letter is to note that elementary 
deductions of the properties of neutron-star fault 
planes show that brittle fracture is not possible. 

2. Neutron-star fault planes 

The relation between stability and stratification 
for neutron star matter (Reisenegger & Goldreich 
1992; see also Jones 2002) constrains the move- 
ment of matter, bounded by any fault plane which 



may be formed, to an almost exactly spherical 
equipotential surface. Shear strain (for in-plane 
shear) is shown in Figure 1 by the displacement 
of a series of constant surfaces in Lagrangian co- 
ordinates which intersect a local element of the 
fault plane xz. In the neutron star case, it might 
be caused by a changing Maxwell stress reflect- 
ing evolution of the internal magnetic flux den- 
sity distribution. The components atj of the stress 
tensor act on the surfaces of the volume elements 
shown on opposite sides of the fault-plane. (We 
assume an isotropic elastic medium in which, be- 
fore fracture, stress propagates perfectly across the 
plane.) In the brittle fracture of a terrestrial earth- 
quake, the stress falls at the instant of failure from 
/ = o~ xy to a much smaller value (zero in the ideal 
case) . Stress energy is largely converted to kinetic 
energy so that strain relaxation Ae occurs with 
an acceleration v such that shear waves are gener- 
ated efficiently (see Kostrov & Das 1988 concern- 
ing the definition of a tectonic earthquake source) . 
The mechanical properties of terrestrial matter, 
as a function of depth, have been tabulated, for 
example, in Kaye & Laby (1986). At a typical 
earthquake focus depth of 15 km, the pressure 
P w 10 _2 /i, where \x is the shear modulus. Crack 
propagation is possible because P is not so large 
that it inhibits void formation behind the tip. At 
greater depths, the increases in T/T m , where T m is 
the melting temperature, and in P/fi both inhibit 
the void formation which is necessary for crack 
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propagation. Thus a transition to plastic stress- 
response occurs at approximately 20 km depth be- 
cause dislocation glide is not pressure-inhibited 
(see Scholz 1990). Tectonic earthquake sources 
are not found in the plastic region. The deep- 
focus sources, which occur in subduction zones at 
depths up to 500 km, are probably intermittently 
running polymorphic phase transitions, following 
the original idea of P. W. Bridgman (1945). 

The contrast with neutron-star matter is ex- 
treme because the latter is not absolutely sta- 
ble, being in equilibrium only at finite pressure. 
Those isotropic components of the stress tensor 
derived from the electrons and from Coulomb in- 
teraction (the Coulomb-electron partial pressure 
Pee which excludes the neutron partial pressure) 
are one or two orders of magnitude larger than 
the shear modulus and this largely determines the 
structure of defects such as monovacancies (Jones 
1999). The way in which the nearest neighbours 
to a monovacancy site relax by displacement con- 
siderably reduces the monovacancy formation en- 
thalpy, for example, to 13 Mev at a matter density 
of 8x 10 13 g cm" 3 in the neutron-drip region. Such 
relaxation is not possible for the nuclei surround- 
ing a void whose linear dimension is of the order 
of, or greater than, the relativistic electron screen- 
ing length. In any solid for which the inequality 
Pee ^> H is satisfied, it is anticipated that a void, 
even if formed, would be unstable against dissocia- 
tion to monovacancies with a lifetime, in neutron- 
star matter, perhaps of the order of 10 -21 s. Thus 
we can assert that nuclei at a fault plane may 
have different, probably less, small-distance order 
than exists elsewhere but almost the same nearest- 
neighbour separations. Reference to the theory of 
crack stability (Anderson 1995, see also Landau & 
Lifshitz 1970) shows that, in brittle fracture, the 
important properties of the stress distribution of 
an unstable crack are the zero on the fault plane 
and the singularity at its tip (in the ideal case) 
which enable it to propagate into regions of much 
lower stress. These do not exist in neutron stars 
owing to the impossibility of forming a sufficiently 
long-lived void. 

Previous formation enthalpy calculations (Jones 
2001) have shown that an amorphous heteroge- 
neous solid phase is formed in the crust and is 
likely to persist as the star cools. The mean square 
deviation (AZ) 2 from the average nuclear charge 



is large compared with unity. The system is anal- 
ogous with a multi-component amorphous alloy 
in having some local order, which vanishes with 
increasing length scales, but differs in having no 
absolute stability. Shear stress acting on such a 
system produces faults in the form of localized 
shear bands. These are thin layers of inhomo- 
geneuos plastic flow, formed sequentially, which 
permit local strain relaxation through through 
dissipation and the transfer of stress energy to 
neighbouring regions. A shear band has structure 
similar to a viscous layer being deformed adiabat- 
ically between parallel plates (see, for example, 
Xing ct al 1999, who also give a typical stress- 
strain relation). The process is not analogous with 
brittle fracture because most of the local stress en- 
ergy, instead of being converted to kinetic energy, 
is either dissipated or transferred to neighbour- 
ing regions. The rotation of a circular cylinder of 
solid crust under Maxwell stress, by plastic flow 
on the cylindrical fault surface delimiting it, is 
an elementary example satisfying the stratifica- 
tion and stability constraint (sec also Thompson 
ct al 2000). The plastic flow threshold is f a and 
the Maxwell stress is assumed to be / exp(— s/s ), 
decreasing with slip distance s, where s Q and / 
are constants. Suppose that an unspecified stress- 
transfer incident results in Af = f — f a > at 
some instant. The equation of motion is 

-pas = Af - f- (1) 

4 s 

where p is the matter density and a the cylin- 
der radius. The velocity at time t during slip 
through a distance 2s p , where s p — s Af/f, 
is v = s = ujs p sin((2)t). The angular velocity 
uj = y/4f/pas is small compared with the values 
2q4-5 ra( j g -i £ Qr wn j cn there is efficient coupling 
of shear and Alfvcn waves. ( It is Co = 6.3 x 10 2 
rad s" 1 for the maximum shear stress 10 27 dyne 
cm" 2 , p = 10 14 g cm" 3 , and s a = a = 10 4 cm.) 
In principle, adiabatic heating on the fault surface 
could lead to a localized solid-liquid phase transi- 
tion, the reduction in viscosity by many orders of 
magnitude giving slip velocities much larger than 
Cjs p . An order of magnitude estimate of the con- 
dition that slip and adiabatic heating should in- 
crease the temperature to T m 

on the fault surface 
is given by f s p ~ CT m yJ(n/ujC), where k is the 
thermal conductivity and C the specific heat. At a 
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typical density of 8 x 10 13 g cm" 3 , T m = 5.8 x 10 9 
K for the lattice parameters given by Negele & 
Vautherin (1973). The specific heat is almost en- 
tirely that of normal neutrons; C — 4.9 x 10 20 
erg cm _3 K _1 for an effective mass of 0.8m„ (see 
Pines & Nozieres 1966). The thermal conductiv- 
ity k = 1 x 10 20 erg cm _1 s _1 K _1 has been esti- 
mated from the electron conductivity results given 
by Gnedin ct al (2001). For f Q w 10 27 dyne cm" 2 , 
this condition shows that a local slip distance ex- 
ceeding s p w 50 cm would be required to increase 
the fault surface temperature to T m . Thus a slip 
of this order of magnitude could produce melting 
and subsequent slip velocities much larger than 
uosp. But the assumption of a single shear band in 
the form of a cylindrical fault surface is unique and 
does not appear realistic. The general case is of 
more complex movements, with considerable de- 
formation, in the form of plastic flow where stress 
transfer produces many distinct shear bands with 
small slip distances not satisfying the melting con- 
dition. In such flow, the velocity v w c s Ae and ac- 
celeration v « lov necessary for efficient shear wave 
generation at angular velocity ui are not reached. 

3. Conclusions 

Our conclusion is that the strain relaxation con- 
ditions necessary for the generation of shear waves 
are not present in neutron stars. This need not 
have a great impact on theories of the pulsar glitch 
phenomena. In the theory of Ruderman, Zhu & 
Chen (1998), the movements caused by Maxwell 
stress and by spin-down of the neutron super- 
fluid in the core can be many orders of magnitude 
slower than those associated with brittle fracture 
and yet remain consistent with the observed up- 
per limit on the spin-up time for Vela glitches. But 
the hypothesis that a plastic-brittle transition, on 
cooling to internal temperature T w 10 _1 T m at 
ages between 10 3 and 10 4 yr, is responsible for the 
absence of large glitches in very young pulsars is 
not consistent with the properties of neutron-star 
matter fault-planes deduced here. These may also 
be useful in determining the mechanism for en- 
ergy release in AXP and SGR sources (Thompson 
& Duncan 1996, Thompson et al 2000). These au- 
thors assume that a plastic-brittle transition oc- 
curs at internal magnetic flux densities near to 

= y/^njl, and that brittle fractures at B < 
arc the source of Alfven wave generation (Blacs 



et al 1989) leading to both quiescent and burst 
emission. But the mode of fracture is determined 
by the nature of the solid, and its occurrence de- 
pends not on the absolute size of Maxwell stress 
components but on the extent to which they de- 
viate from values that would give an equilibrium 
in a completely liquid star. Our conclusion, that 
brittle fracture is not possible at any temperature, 
indicates that other mechanisms (see, for example, 
Thompson, Lyutikov & Kulkarni 2002; Mereghetti 
et al 2002) need to be considered in more detail. 
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Fig. 1. — For in-plane shear, with displacement 
shown by a series of constant surfaces in La- 
grangian coordinates, the components of the stress 
tensor acting on volume elements on opposite sides 
of the y = fault-plane are / = a xy before fail- 
ure. In ideal brittle failure, / = immediately 
afterwards and strain relaxation occurs with ac- 
celeration v such that stress energy is efficiently 
converted to shear waves. For neutron-star mat- 
ter, / can exhibit no sudden decrease. 
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